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The critical behavior of Ising model on a one-dimensional network, which has long-range con- 
nections at distances I > 1 with the probability O(Z) ~ l~ m , is studied by using Monte Carlo 
simulations. Through analyzing the Ising model on networks with different m values, this paper 
discusses the impact of the global correlation, which decays with the increase of m, on the phase 
transition of the Ising model. Adding the analysis of the finite-size scaling of the order parameter 
[(M)], it is observed that in the whole range of < m < 2, a finite-temperature transition exists, 
and the critical exponents show consistence with mean-field values, which indicates a mean-field 
nature of the phase transition. 
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I. INTRODUCTION 

Since first proposed by Watts and Strogatz [l|] in 1998, 
the Watts-Strogatz (WS) small world model has been 
widely and deeply studied. In WS model, vertices are 
placed on a ring with each vertex having a finite num- 
ber of 2k nearest regular connections initially. The con- 
nections are then rewired with a probability p r to form 
long-range connections or short paths. By varying the 
single parameter p r from to 1, WS model displays phase 
transition from a regular network to a small world, which 
ends up with a random network at p r — 1. Considering 
that the rewiring of connections may cause isolated ver- 
tices, the addition-type WS model - Newman and Watts 
(NW) model 0, Q was proposed later, in which long- 
range connections are added with a probability p while 
keeping the initial regular connections unchanged. NW 
model and WS model are nearly equivalent, especially 
when the network size N — > oo @, H. Marcelo Kuper- 
man and Guillermo Aberamson generalized WS model 
in Q. In their model, besides connecting with the two 
nearest neighboring vertices, with a probability p each of 
the vertices in the model chooses a non-neighboring ver- 
tex to construct a long-range connection between them. 
The choice of the non-neighboring vertex is governed by 
the distance-dependent probability distribution: 
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When m = 0, the NW model is restored. Through 
varying the value of to, the global correlation of the net- 
work is influenced. By analyzing the topological charac- 
teristic of the network, it was observed that the network 
also shows a small world nature, and the network con- 
verges to a regular lattice as to — > oo 0, H, H, 0] • 

To analyze the behavior of a small world network from 
the point of view of statistical physics, researchers focus 
on its long-range connections or short-paths, which may 



lead to global coherence. In networks without any short- 
path, the transmission of information has to pass a long 
distance ~ O(N) (N is the network size), so global coher- 
ence is difficult to reach. With the increase of the number 
of short-paths, or long-range connections, collective be- 
havior will emerge in the system. This global coherence 
and the ubiquitous small world phenomenon in real-life 
make the study of thermodynamic phase transition on 
networks with small world property significantly popular 
H ©, E3, El E3, El, in which Ising model is one of the 
most fascinating points. 

As a comparatively simple but very important model 
of statistical physics, Ising model perfectly shows order- 
disorder transition of the system. Ising model on one- 
dimensional regular lattice does not show phase transi- 
tion at any finite-temperature and the phase transition 
on multi-dimensional regular lattice is not of mean-field 
nature. That the Ising model undergoes phase transi- 
tions in the addition-type WS model has been studied in 

nans 

, [13, HH • It has also been shown that in the small 
world phase the addition-type WS model has a mean- 
field nature [TJ E3, EH- The critical behavior of the 
Ising model on one-dimensional network with distance- 
dependent connections given by Eq.(l) is expected to re- 
flect the nature of the network at different values of m 
indirectly But it remains controversial about the 

mean-field nature of the phase transition in the range of 
1 < to < 2 0,0,113. 

To detect whether a network has a small world be- 
havior, one can calculate the average shortest-path 
length and the average clustering coefficient. Sen and 
Charkrabarti in ((J and Moukarzel and de Menezes in Q 
obtained the contradictory results by analyzing the topo- 
logical structure of networks with connection probability 
at distance I given by Eq.(l). In Ref. [6(, it was found 
that the average shortest-path length behaves as In N on 
rings of size N for all m < 2, hence it was concluded 
that small world behavior occurs for < to < 2 while in 
Ref. 0] , it was argued that small world behavior occurs 
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only for m < 1. In the section 1 < m < 2, according 
to Ref. Q, the average shortest-path length scales as 
N s with the value of the exponent 5 (0 < S < |) de- 
pending on ui. Ref. [l2| studied the critical behavior of 
Ising model on such a one-dimensional network, since if 
the network's behavior is that of small world it should 
be reflected in the critical exponents of the Ising model, 
which will assume mean-field values. Their result indi- 
cates that there is a mean-field behavior for m < 1 and 
a finite-dimensional behavior for 1 < m < 2, which is in 
accordance with the conclusion in Ref. Q- 

In this paper, we re-examine the critical behavior of 
Ising model on a one-dimensional network with distance- 
dependent connections in Ref. [13] by using Monte Carlo 
simulations. But we enlarge the underlying network sizes, 
and also take into account the finite-size scaling analy- 
sis of the order parameter [(M)] by using the plot of 
[(M^/V 1 / 4 versus T, which can reflect the mean-field na- 
ture indirectly [8j. Here, [• ■ • ] and (■ • • } denote the ther- 
mal average taken over Monte Carlo steps for equilibrium 
at each temperature, and over different network realiza- 
tions, respectively. Thus more reliable results can be pro- 
duced from our simulations, since statistical fluctuations 
are greatly suppressed in these large-scale networks, and 
the finite-size effects are much more discernible in the 
Binder's cumulant Un [H, H3] and the specific heat C v 
than in the order parameter [(M)\ Q . The results show 
that the model has a mean-field nature in the whole range 
of < m < 2. 

In the next section, we will introduce the model we use 
and the Monte Carlo simulation. In section III, the re- 
sults and analysis are given. Section IV is the conclusion. 



II. ISING MODEL AND MONTE CARLO 
SIMULATIONS 




FIG. 1: Examples of generated networks with TV = 101 under 
the same p — 1 and different m values. 



regular network with 2k = 4 nearest neighbors. Fig.l 
shows that with the increase of m, long-range connections 
are more and more difficult to form between vertices with 
long distances. 



The underlying network of the one-dimensional Ising 
model used in this paper is based on the generalized WS 
model [H . The addition of connections considers the dis- 
tance between vertices. That means the vertex j to which 
a connection is attached is not randomly selected, but ac- 
cording to a distribution depending on the distance I from 
i to j: 0(1) ~ l~ m . We start with a one-dimensional reg- 
ular lattice with N = 2K + 1 vertices and 2k = 2 nearest 
neighbors. From i = 0, we select the i-th vertex and 
generate a random number tp E [0,1) from a uniform dis- 
tribution. If <p < p, we select one vertex from the clock- 
wise K vertices to attach to according to the probability 
distribution 0(1) ^ l~ m . Self-connections and multiple 
connections are prohibited. We repeat this process until 
all the vertices are selected, which results in a network 
with N(l + p) connections on average. 

In order to investigate the influence of m on the critical 
behavior of Ising model, in this paper we only consider 
the case p = 1 for simplicity. Our model with m = 
is equivalent to the model used by Andrzej Pekalski Q. 
As m — ► oo, our model restores to the one-dimensional 
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FIG. 2: The relationship of the average shortest-path length 
L with the size of the network N under different m values 
(the case m = 0.5 is almost identical with m — 0). All the 
plots have been averaged over 50 realizations. 

Fig. 2 shows that the increase of m results in the in- 
crease of the average shortest-path length L, but the 
small world property retains within the entire range 
< m < 2 as L oc InN. Fig. 3(a) shows that with the 
increase of m value, especially after m > 1, the network 
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clustering coefficient becomes independent of system size, 
which is the characteristic of regular network. Further- 
more, it can be seen from Fig. 3(b) that, for a certain 
network size, with the increase of m, the clustering coef- 
ficient converges to C max = 1/2, which is identical to the 
clustering coefficient of one-dimensional regular network 
with 2k = 4 nearest neighbors. 



1=0.0 




temperature, and (• ■ • ) denotes the average over different 
network realizations taken over 30 — 80 configurations. 

Besides these three physical quantities, the finite-size 
scaling analysis of the order parameter [(M)], which ex- 
hibits the critical behavior [{M)\ ~ {Tc-T) 13 in the ther- 
modynamic limit, is used. In a finite-sized system, the 
order parameter scales as Q: 



\{M)\ = N-P' v g{{T - TJN 1 ^), 



(4) 



where g(x) is an appropriate scaling function. Eq.(4) 
leads to a unique crossing point at T c in the plot of 
\(M)]N f3 / v versus T, so Beom et al Q suggested ana- 
lyzing the finite-size scaling of the order parameter by 
[(M)]^ 1 / 4 for large m to overcome the finite-size effects, 
which are more prominent in other thermodynamic quan- 
tities. 



FIG. 3: (a) Network clustering coefficient C versus network 
size N for different m values, (b) Network clustering coef- 
ficient C versus m value for a network with TV = 100001 
vertices. The shown C values in these two plots have been 
averaged for 50 different network realizations. 



III. RESULTS AND ANALYSIS 



The Ising model is described by the Hamiltonian: 

(i,3) 



(2) 



where J > is the coupling constant between vertex i 
and j if they are connected, (i, j) is the collection of all 
the connections in the network. 

Our Monte Carlo simulation of the Ising model starts 
on a periodic one-dimensional ring. At the beginning, 
all the spins on the ring are given the same value +1, 
then a random vertex is selected to flip according to the 
rule of the Metropolis algorithm 0, [l8| and Glauber dy- 
namics [ll| to simulate the evolution of Ising model un- 
der different temperature. And finite-size scaling is used 
to investigate the paramagnetic-ferromagnetic transition 
temperature and the critical exponents. During the sim- 
ulations, we mainly use the random number generator in 
[l9| . and we also use the drand48() function in standard 
C library as a comparison. It is found that there is no 
notable influence on the simulation results. In the simu- 
lations we measure the Binder 'cumulant Un, the suscep- 
tibility x and the specific heat C v [H, [H 0, HH, HI : 



U N = 1- 



[(M 4 )] 
3[(Af 2 )] 2 ' 
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(3) 



m 2 ) - (h) 2 ] 



with M 



T 2 N 

Here [■ • • ] denotes the thermal av- 



erage taken over 2.5 x 10 4 Monte Carlo steps after dis- 
carding the initial 2.5 x 10 4 ones for equilibrium at each 
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FIG. 4: For m = 0, (a) Binder's cumulant Un has a unique 
crossing point at T c = 3.10(5)(in units of J/fcs). (b) The 
critical exponent V = 2.08 is obtained from the linear fit of 
InAf/jv vs. In A. (c) Specific heat C v has a unique cross- 
ing point at T c = 3.10(7), which indicates a = 0. (d) The 
expansion of C„ near T c obtains V = 2.04. (e) Finite-size 
scaling of the susceptibility again determines T c — 3.10(5) 
with the critical exponent 7 = 1. (f) Finite-size scaling of 
the order parameter [{M)] leads to a unique crossing point at 



T c = 3.10(1), and yields /3 = 



2.0 (see Eq.(4) |g). 



As shown in Fig. 4, for m = 0, the finite-size scaling 
of the Binder's cumulant Un, the specific heat C v , the 
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susceptibility \ an d the order parameter [(M)] all re- 
veal unanimously T c = 3.10(5) (in units of J /ks), which 
confirms the presence of a finite-temperature transition. 
And the obtained critical exponents indicate the mean- 
field nature of the transition. 

Till now, the most controversial property of the phase 
transition of one-dimensional Ising model with distance- 
dependent connections given by Eq.(l) is in the range of 
1 < m < 2 H 0, [T§], so we will focus on the critical 
behavior of Ising model for 1 < m < 2. 
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FIG. 5: For m = 1.0, (a) Finite-size scaling of the Binder's 
cumulant Un has a unique crossing point at T c = 3.05(2). 
(b) Finite-size scaling of the specific heat C v determines T c = 
3.05(3), which indicates a = 0. (c) The expansion of C v 
near T c obtains V — 2.05. (d) Finite size-scaing of the order 
parameter [{M)] leads to a unique crossing point at T c — 
3.05(3), and yields /3 = §, P = 2 (see Eq.(4) ||). 

Fig. 5 presents the finite-size scaling of the Binder's cu- 
mulant U n , the specific heat C v and the order param- 
eter [(M)] at m — 1.0. The measured quantities re- 
main a unique crossing point and reveal unanimously 
T c = 3.05(3), which confirms the presence of a finite- 
temperature phase transition. The obtained critical ex- 
ponents assume the mean-field values as well. 

Fig. 6 is the finite-size scaling of the Binder's cumulant 
Un, the specific heat C v and the order parameter [(M)] 
at to = 1.5. The curves of Binder's cumulant Un do not 
intersect to one unique crossing point very well. But the 
trends of the evolution indicate that it will emerge when 
the underlying network size is large enough, the finite- 
size scaling of the order parameter \(M)] under different 
network sizes behaves better in intersecting to one unique 
point, revealing unanimously T c ~ 2.70(2). The obtained 
critical exponents values indicate the presence of a mean- 
field finite-temperature phase transition. 

With the increase of to, the algebra decreasing of 
0(Z) ~ l~ m will will have a cut-off due to the finite net- 
work size. This makes the normalization of the distri- 
bution function depend on the size of the network, and 
makes it inevitable to study systems of very large sizes 




FIG. 6: For m = 1.5, (a) Finite-size scaling of the Binder's 
cumulant Un has a unique crossing point at T c — 2.70(2). (b) 
The expansion of Un near T c obtains V — 2.03. (c)Finite-size 
scaling of the specific heat C v determines T c — 2.70(4), which 
indicates a — 0. (d) Finite size-scaing of the order parameter 
[(M)} leads to a unique crossing point at T c — 2.70(2), and 
yields (3 = ±, v = 2 (see Eq.(4) @|). 
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FIG. 7: Finite-size scaling of the order parameter [(A/)] gives 
T c = 2.58(5) at m = 1.6 (a) and T c = 2.28(2) at m = 1.8 (b). 



so as to obtain correct scaling behavior [8|, [22j . By an- 
alyzing the plot of [(M)JTV 1 / 4 versus T, we get further 
progress on the transition temperature of systems with 
large scale, as shown in Fig. 7 as examples. This also 
indicates the fact that the phase transition still holds a 
mean- field nature for to > 1. The phase transition tem- 
peratures here are higher than those found by H. Hong et 
al in Ref. [20| , which may be due to the presence of more 
long-range connections in our model which contribute to 
different phase transition properties. 

In Fig. 8, the decrease of T c in the range of [0,1) is 
not significant, but this decreasing becomes obvious af- 
ter to > 1 . This brings difficulties to the phase transition 
analysis of Ising model due to the following two factors: 
a) With the increase of to, the appearance of a convinc- 
ing unique crossing point needs larger network sizes, as 
shown in Fig. 6 and Fig. 7; b) The transition tempera- 
ture decreases rapidly after to > 1 (Fig. 8), so we need to 
do Monte Carlo simulation at lower temperature, which 
means lower probability to flip. What's more, the value of 
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FIG. 8: The evolution of the transition temperature with the 
increase of m. 

the order parameter [(M)] has a strong correlation with 
its value at the previous step. Thus the system needs 
more Monte Carlo steps to reach equilibrium, and we 
need more iterations to get the thermodynamic average. 
Thus the simulations will be more time-consuming. 

IV. CONCLUSION 

In this paper, we reconsider the phase transition of 
Ising model on a one-dimensional network with distance- 
dependent connections given by 0(7) ^ l~ m , by adding 
the finite-size scaling analysis of the order parameter 



\(M)] to overcome the difficulties at high to values. 
Ref. [12| based their model on networks with hundreds 
of vertices and concluded that the phase transition has 
a mean- field nature only in the range of < to < 1. 
Our results, based on larger network sizes and further 
finite-size scaling analysis of the order parameter [(M)] 
as well as the Biner's cumulant £/jv, the specific heat C v 
and the susceptibility x, show that the phase transition 
has a mean- field nature in the whole range of < to < 2. 

Because of the constraint of the long simulation time, 
we can not apply the analysis to networks with even 
larger sizes at present. Thus we can not make sure 
whether there is a critical value of to at or beyond to = 2, 
above which there will be no mean-field nature or finite- 
temperature phase transition since the transition point 
to = 2 was suggested by the study of the network topo- 
logical structure [1, 0, GjJ • With the finding of a proper 
solution to the long simulation time, more interesting re- 
sults are supposed to come out. 
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